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Abstract 

We study, analytically and numerically the stationary states in the system of two linearly coupled 
nonlinear Schrodinger equations in two spatial dimensions, with the nonlinear interaction coeffi- 
cients of opposite signs. This system is the two-dimensional analog of the coupled-mode equations 
for a condensate in the double- well trap [Physical Review A 69, 033609 (2004)]. In contrast to the 
one-dimensional case, where the bifurcation from zero leads to stable bright solitons, in two spatial 
dimensions this bifurcation results in the appearance of unstable soliton solutions (the Townes- 
type solitons). With the use of a parabolic potential the ground state of the system is stabilized. 
It corresponds to strongly coupled condensates and is stable with respect to collapse. This is in 
sharp contrast to the one-dimensional case, where the ground state corresponds to weakly coupled 
condensates and is unstable. Moreover, the total number of atoms of the stable ground state can 
be much higher than the collapse threshold for a single two-dimensional condensate with a negative 
s-wave scattering length. 
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I. INTRODUCTION 



Bose-Einstein condensates (BECs) in trapped dilute gases exhibit interesting interplay 
between quantum coherence and nonlinearity since, at zero temperature, the quantum gas 
is described by the mean-field theory based on the nonlinear Schrodinger equation with an 
external potential - the Gross-Pitaevskii (GP) equation |l| for the order parameter. The 
macroscopic quantum coherence of BEC, first experimentally demonstrated in Refs. 
was subsequently explained theoretically |4] with the use of the GP equation. 

Nonlinear phenomena in BEC bear similarity with nonlinear optics. Similar to optics, 
where the bright and dark solitons are supported respectively by the focusing and defocusing 
nonlinearities, in BECs the s-wave scattering length is the determining factor. Dark soli- 
tons are routinely observed in the in the quasi-one-dimensional condensates with repulsive 
tactions On the other hand, the attractive W ln _, BEC plates 

in the form of the bright solitons P. llo|. 

The s-wave scattering length can be modified by application of magnetic field near the 
Feshbach resonance For a BEC constrained to lower spatial dimensions, the Feshbach 



resonance still proves to be sharp, for instance, in the two-dimensional condensate 



121. 



The feasibility of control over the scattering length in BEC by the optical means was also 



proposed |l3( (see also Refs. |l4 Il5j|). 



Control over the scattering length in a part of the condensate can be realized in the double- 
well trap with far-separated wells. Condensates in the double-well potential are currently 
routinely created and studied in the experiments (see, for instance, Refs. 16, 17, 18^. 

The double-well trap is created in one spatial dimension, the other two dimensions thus 
allow for two geometrically distinct cases which correspond to the one- and two-dimensional 
BECs depending on the trap asymmetry. A combination of the double-well trap with the 
control over the scattering length allows one to observe two tunnel-coupled BECs with the 
opposite interactions (i.e., one attractive and the other repulsive). In on e sp atial dimension 
such a setup leads to appearance of the unusual stable bright solitons [3 Q which have 
almost all atoms contained in the repulsive condensate. The one-dimensional ground state, 
lowever, corresponds to weakly coupled condensates and is unstable with respect to collapse 
3]. In this paper we study the two-dimensional case. We find that the two-dimensional 
solitons (the Townes type solitons) are always unstable. However, with the help of a confining 
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potential, the ground state can be stabilized. In contrast to the one-dimensional case, in the 
two- dimensions the ground state of the system corresponds to strongly-coupled condensates 
and there is an energy barrier for collapse (see also Ref. |21|.) 

The applicability of the GP-based mean-field theory is limited, but it always applies 
to the description of stable stationary states. This is due to the similarity between the 
Bogoliubov-de Gennes equations, describing the surrounding cloud of hot atoms, and the 



221. Thus, a 



equations describing evolution of a linear perturbation of the order parameter 
stable mean-field stationary state is also stable in the full quantum approach, its life time 
is equal to that of the condensate. 

Theoretical investigations of BEC in the double- well trap go back to the prediction of 



the anomalous 
condensate 24 



25 



Experimental advances 



osephson oscillations 
23. 

Q i 



and the macroscopic quantum self-trapping of the 



16 



17 



in the production and manipulation of the conden- 
sate in the double-well trap make the implementation of the tunnel-coupled repulsive and 
attractive condensates feasible. Recently a direct observation of the quantum tunneling and 
nonlinear self-trapping of a BEC in the double- well trap was reported [lsj ]. 

The collapse instability in an attractive BEC is loaded into the double-well trap was also 
recently investigated. In the quasi-one dimensional case, the critical number for collapse 
;urns out to be larger than the same for the corresponding axially-symmetric harmonic trap 
2t] |. The tunnelling-induced collapse of an attractive BEC in a double- well trap can take 
place under the influence of a time-dependent potential 

Control over the scattering length in one part of a condensate can be realized in a double- 
well trap with far separated wells. In this case a simplification of the GP e quat ion is possible, 



which results in the coupled-mode approximation similar to that of Refs. [25j,|26j]. However, 
in contrast to the latter works, the kinetic energy of the condensate is taken into account 
(see secti on ITT1) . The coupled-mode system also has applications in optics (see, for instance, 
Refs. HOB)- 

Since the two-dimensional (2D) nonlinear Schrodinger equation (NLS), 



i<9^ + V 2 ^-#|^| 2 ^ = 0, 



(1) 



is critical, we adopt the point of view based on the analysis of the critical scaling and its 
perturbations. The NLS equation has the following scale invariance: if $?(t, r) is a solution 
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then ty(t,r) = k^(k 2 t,kr) is also a solution. The number of atoms N (or the Z 2 - n orm), 
defined by N = J d ra r|\I/| 2 , is scaled as N = k 2 ~ n N in n spatial dimensions. 

The scale invariance of the 2D NLS equation allows for a family of solutions with the 
same number of particles. One may call the 2D scale invariance the "critical scaling". 
The critical scaling leads to important consequences (see, for instance, Ref. 32 1). As the 
number of particles is constant for the whole family of solutions, the Vakhitov-Kolokolov 
(VK) criterion 33] applied to the Townes soliton gives marginal stability (or instability), 
since dN/dfi = 0, where \i is the chemical potential for a particular solution of the family 
(i.e., —ft is the frequency). 

This explains why addition of an external confining potential allows for stable localized 
solutions. Indeed, the external potential breaks the scale invariance and the number of 
particle degeneracy is broken too: some solutions of the former Townes soliton family have 
number of particles below the collapse threshold (moreover, thanks to the so-called lens 

n 

transformation 34] , the collapse threshold does not depend on the strength of the potential 
if the latter is parabolic). 

Besides addition of an external potential, there are other ways to break the scale in- 
variance. The coupled-mode system, i.e. the system of two linearly coupled NLS equations, 
provides another way. The coupling coefficient is proportional to the tunnelling rate through 
the central barrier. 

Having understood the relation between the broken scale invariance and stability against 
the collapse, one may wonder if the linear coupling of two 2D NLS equations allows for the 
existence of stable two-dimensional (i.e. Townes-type) solitons. Indeed, if the critical scaling 
is broken and the number of atoms (or the number of particles, generally) depends on the 
chemical potential then there could be self-localized stationary solutions (which have right 
to be called solitons) with dN/dfi < 0, i.e. satisfying the VK criterion, with the hope that 
they are stable. 

In the following, using the singular perturbation theory supplemented by numerical sim- 
ulations, we argue that the two-dimensional soliton solutions to the coupled-mode system 
are always unstable with respect to collapse. The reason is that, in contrast to the ID 
case , the bifurcation from zero in the 2D coupled- mode system is always discontinuous: 
the number of atoms of the two-dimensional soliton solution with vanishing amplitude does 
not vanish. However, addition of an external confining potential (a potential transverse to 
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the double-well trap) restores the continuity of the bifurcation from zero and leads to the 
appearance of stable solutions. Some of them have the number of atoms much larger than 
the collapse threshold in a single 2D NLS equation. Moreover, such a solution is the ground 
state of the system, since it has the lowest possible energy for a fixed number of atoms. 
In contrast to the one-dimensional case, the ground state in the two-dimensional system is 
secured from collapse by an energy barrier. 

The paper is organized as follows. In section [H] we derive the coupled-mode system 
from the Gross-Pit aevskii equation for a condensate in an asymmetric double- well trap. In 
section IIHI we consider the stability properties of the axially symmetric stationary states. 
Section HVI is devoted to the study of the limiting case solutions: the bifurcation from zero 
in the coupled-mode system with and without the parabolic potential, sections II V Al and 
IIVBI respectively, and the asymptotic solution corresponding to large negative values of 
the chemical potential, section HV CI The reader not interested in details of the bifurcation 
analysis may go directly to section where we summarize the results and discuss the 
numerical solution of the coupled-mode system. In the numerics we have used the Fourier 
spectral collocation method and looked for the stationary solutions using the numerical 
schemes of Refs. JsjJ and j^]. Section IS contains some concluding remarks. 



II. DERIVATION OF THE COUPLED-MODE SYSTEM 

The coupled-mode system follows from the GP equation under two conditions. First, 
the wells of the double-well trap must be far separated. Second, the number of BEC atoms 
must be below a certain threshold (equation (fTTj) below) with the result that the motion of 
the condensate in the spatial dimension of the double-well trap is equivalent to that of a 
quantum particle. The stationary states considered below satisfy this condition. It follows 
that BEC atoms occupy only the degenerate energy levels of the double-well trap. 

The tunnelling coefficient, usually defined through an integral over the overlap of the 
wave functions (see, for instance, Ref. SslE^)' can be, in fact, explicitly given in terms of 
the parameters of the double- well trap (equation (jHJ)). This is due to the above conditions 
and the simple fact that the wave functions of the localized basis (given by equations (jSJ) 
and (jEJ)) are uniquely defined by the trap. 

The two-dimensional coupled-mode system describes the so-called pancake condensates. 
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Although the ID and 2D coupled-mode systems are similar, there is a difference in the order 
of magnitude of the respective parameters which is explained below. 

The Gross-Pitaevskii equation for the order parameter ty(r,t) of a BEC in a double- well 
trap given by a parabolic potential with a Gaussian barrier reads 

ihdt*= {-^V 2 + V ext (r)+g\V\ 2 ^, (2) 



where 

V^^^lf^ + ^ + Vsexpl- 1 -^^}. V U >0. (3) 



Here the parameters zq and a give the position and width of the barrier, respectively. The 
pancake geometry corresponds to the condition 7 = ou z /lux ^ 1- I n terms of the potential 
from equation (JSJ) the wells of the double well trap are far separated if a ~ £ z with £ z = J 
being the oscillator length in the ^-direction. 

The motion in the double- well trap and the transverse dynamics of BEC can be factorized 
if the number of BEC atoms is not large (see equation (jlljl ). i.e. if the oscillator length i z is 
much smaller than the characteristic length of nonlinearity £ n \. In the case of a double- well 
with far separated wells, the first two energy levels are quasi-degenerate: 

E 1 -E <^E 2 - El (4) 

Being interested in the ground state of the system (occupied by the condensate) we consider 
only the degenerate subspace. This approximation neglects the uncondensed atomic cloud 
which can be discarded for temperatures below the condensation threshold 

The basis in the degenerate subspace for expansion of the order parameter is dictated by 
the necessity to simplify the nonlinear term, which is not small for the transverse degrees of 
freedom, i.e. in the pancake plane. Evidently, one should select such a basis, say ip u {z) and 
ip v (z), where each basis function is localized in just one of the wells. The localized basis is 
given by a rotation of the wave functions for the ground state and the first excited state: 

V 1 + X V 1 + X 1 

Obviously, the new wave functions are orthogonal and normalized. The parameter x is 
selected by the quotient of the absolute values of the eigenfunctions ipo(z) and ij)\{z) at the 
minima of the double-well potential z_ and z + (say, z_ < z + ) [l^. We set 



(the positions of the extremals of the wave functions slightly deviate from the minima of the 
trap; these deviations we neglect). For this choice of x the wave functions ip u (z) and ip v (z) 
defined by equation (jSJ are localized in the left and right well, respectively (see also figures 
1 and 2 in Ref. H)- 

The Hamiltonian for a quantum particle in the double-well potential, projected on the 
degenerate subspace, i.e. H z = Eo\ip ) (ip \ + Ei\ipi)(ipi\, in the new basis becomes 

H z = E(\ip u )(ip u \ + \ip v )(H) + \ipv)£(ipv\ ~ K(\ip u )(ip v \ + \ip v )(ip u \). (7) 

Here 



We can set E = without loss of generality. Equation (jHJ) gives the tunnelling coefficient K 
and the zero-point energy difference S. 

We approximate the solution of the GP equation (j2J) by a sum of the factorized order 
parameters: 

r ± , z) = $ u (t, r±)ip u (z) + r ± )ip v (z). (9) 

Let us now formulate the condition for the factorization in equation ©. We have neglected 
the nonlinear term in the Gross-Pitaevskii equation as compared to the longitudinal kinetic 
term, that is 

!d|»MI*lW, do) 

(for each of the two wells). The wave functions can be estimated as follows \ip\ 2 ~ l/£ z and 
|$| 2 ~ N/d\, with d± being the transverse radius of the condensate and M the number 
of atoms (in the considered well). Using the expression for the nonlinear coefficient in the 
Gross-Pitaevskii equation g = 47rh 2 a s /m where a s is the atomic scattering length, we 
get the applicability condition for the separation of variables as follows 

8ir\a s \£ z Af 



d 2 , 



« 1. (11) 



Condition (JTTJ) must be satisfied by both condensates in the double-well trap. 

The coupled-mode system is derived by inserting expansion (jH} into equation Q, using 
the orthogonality of the basis functions ip u>v , the formulae 

H z ip u = -Kip v , H z ^) v = £4> v - K^pu, 
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and throwing away the small nonlinear terms involving the cross-products of the localized 
wave functions ip u and ip v (see also the discussion of Ref. [id]). One arrives at the system: 

h 2 

ihd t ®u = -7T- Vl$« + V(r ± )<S> u + g u \<S> u \ 2 <$> u - K$ v , (12a) 
2m 

h 2 

ihd t ®v = -7^Vi$, + V{r ± )$ v + {£ + g v \$ v \ 2 )$ v - K$ u . (12b) 
2m 

2 

Here g u>v = f dzg(z)\i/) u>v \ 4 and V(f±) = ^p-rl 2 (by changing the sign of either ip u or ip v 
one can always set K > 0). 

System (|T2"|) is the basis of our approach. Conditions (JH) and (fTTj) are satisfied by all stable 
stationary states considered below. In the search for two-dimensional soliton solutions we 
will neglect the confining transverse potential V(f±) (i.e. when it is considered as flat on the 
scale of the Townes-like soliton solution). This case will be called below "the coupled- mode 
system without the transverse potential". The atomic interaction in the ^-condensate is 
attractive, g u < 0, while in the f-condensate it is externally modified to repulsive, g v > 0. 

For the numerical and analytical analysis it is convenient to use the dimensionless vari- 
ables defined as follows 



T=^t, ?=* t^J-*-. (13) 
2 £± V mu± 



The order parameters are expressed as 



§u = -—u, <f> v = —v, (14) 

where A = 87r|ai u ^| J dz|-^ u | 4 ^ and is the scattering length in the w-condensate. 
The dimensionless system reads 

idxu = (— Vf + p 2 ) u — \u\ 2 u — kv, (15a) 



id T v = (— V| + p 2 ) v + (e + a\v\ 2 )v - ku. 



(15b) 



Here 



M 

I'.S 



/ dz\i/j v \ 
J dz\ip u \ 



2K 2k E x - E 

hw± l + k 2 hw± 



28 _ 1 - x 2 E 1 - E 

fajj± l + x 2 huj_ 



(16) 
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The number of atoms M in the condensate (the /2-norm) is given as follows 

K,v = J d 2 p1$ u ,„| 2 = AN U , , N u = J d 2 pH 2 , N v = J d 2 p\v\ 2 . (17) 

The quantity N U}V will be referred to as the "number of atoms" for short, since we are 
interested only in the relative shares of the number of atoms in the two condensates and 
the ratio of the total number of atoms to the collapse threshold in a single 2D NLS. The 
transformation coefficient A can be estimated as 5 e , z , , it is of order 10 2 — 10 3 for the current 

87r|a s | ' 

trap sizes in the experiments with BECs. 

The tunnelling coefficient k and the zero-point energy difference e of the dimensionless 
coupled-mode system can take arbitrary values. Indeed, from the definition we have 

2K 2K 28 2£ 

hui± ^ t\uj z ' hw±_ ^ hjj z ' 

i.e. there are two unrelated multipliers, the first, 7, is large and the second is small. 

Finally, we can reformulate the applicability condition (JTTJ) for the coupled-mode system 
in the dimensionless variables: 

N u < 7 r 2 , aN v < 7 r 2 , (18) 

where r u>v = /£± is the dimensionless radius of the condensate. The condition (JTBJ) 
is derived by using the transformation fTT]) with the estimate f dz\ip 

u,v\^ ~ l/^z and that 

7 = ^i/e 

We have verified that condition ()18j) is satisfied (for the pancake trap with 7 > 100) by 
the stable stationary states (figures 3,4,6,7 of section IVT) . 

If more than two energy levels of the double-well trap are significantly occupied by the 
condensate, the conditions for applicability of the coupled-mode system (I12|) are violated. 
In this case, one can use the nonlinear coupled-mode approach which results, however, 
in the nonlinearly coupled NLS equations. 



III. STABILITY OF THE STATIONARY STATES 

BEC in a double-well trap can be unstable with respect to collapse if the atomic inter- 
action is attractive (the w-condensate in the notations of the previous section). By setting 
k = in system (fT3j) we obtain for the u-condensate the focusing 2D NLS equation with an 
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external potential. In the simplest case when the potential is parabolic the collapse threshold 
is given by the /2-norm of the Townes soliton: 

N th = 11.69, (19) 



32, 341. The Townes 



since the collapse threshold is independent of the parabolic potential 
soliton is the solution u = e lT R(p) of the 2D NLS equation, i.e. the function R(p) satisfies 

V 2 R + R 3 - R = 0. (20) 

Here and below the operator V is the gradient with respect to p if is it not explicitly 
indicated otherwise by a subscript. In the following we will need the well-known identity for 
the Townes soliton (see, for instance, l^ ) 

d 2 p R 2 = [ d 2 p(Vi?) 2 = \ f d 2 pR\ (21) 



2 _ 

Let us now consider the problem of stability of the stationary states in the coupled-mode 
system. We are interested only in the axially symmetric stationary states, u = e~ l ^ T U{p) and 
v = e~ lflT V(p), where p = \p\. The stability or instability can be established by considering 
the eigenvalue problem associated with the linearized system. Writing the perturbed solution 
as follows 



u 



e~^ T {U(p) + e- tQT U(p)} , v = e"^ T {V(p) + e~ %nT V(p)} , (22) 





where (W,V) is a small perturbation mode with the frequency Q, one arrives at the following 
linear problem for the eigenfrequency: 

tnl = Ai , (23) 

\Vr J \Vi J \ V/ j \V R j 

with 

Here the scalar operators are defined as follows 

L^ = -(V 2 + f/ 2 + /.)+p 2 , = -(V 2 + 3U 2 + f i)+p 2 , 4 W) = -(V 2 -aV 2 + p-e)+p 2 , 

= -(V 2 -3aV 2 + p-e)+p 2 . (25) 
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First of all, the matrix operator Ao is non-negative for positive stationary solutions, i.e. 
satisfying UV > 0. Indeed, the scalar operators on the main diagonal of A can be cast as 
follows 

L ° ~ U U + U' ~ v v + V 

what can be easily verified by direct calculation. Therefore the scalar product of Ao with 
any vector X = (Xi(p) , X 2 (p)) is non- negative: 

(X\A \X) = J d 2 p [X{L^X X + X*L^X 2 - k(X*X 2 + X*X X )} 



> k I d 2 p 



Xi\l~-X2)/~ 



> 0. 



Here we have used the positivity of the operators — h VU 2 Vjj and — -^VV^V^. The operator 
Ao has one zero mode given by the stationary point itself: Z = (U, V), AqZ = 0. 

Non-negativity of A is an essential property for the following, therefore we will concen- 
trate on the solutions satisfying UV > 0, which can be termed "positive" solutions, while the 
ones, satisfying UV < 0, are discarded from the consideration below. The latter solutions 
bifurcate from zero at a higher energy then the positive ones (see sections ITV^l and EB} • 

For a positive solution, the lowest eigenfrequency of the linear stability problem can be 
found by minimizing the quotient 

tf=min < A ;i A f> (2 6 ) 

in the space orthogonal to the zero mode of Ao: (Z\X) = (here (X\Y) = J d 2 p (X*Yi + 
^2X3)). Equation follows from the eigenvalue problem rewritten as A AiX = Q 2 X 
withX = (U R ,V R ). 

The imaginary eigenfrequencies Q, which mean instability, appear due to negative eigen- 
values of the operator Aj_. For the coupled-mode system without the transverse potential 
there is at least one zero eigenvalue of Ai due to the translational invariance, with the eigen- 
function (in the vector form) being given as Z\ = (VU, W). For the positive stationary 
solutions, if there is only one negative eigenvalue, the VK stability criterion ^ < applies, 
which can be established by a simple repetition of the arguments presented, for instance, in 
Ref. (22I]. The limit on the number of negative eigenvalues is related to the fact that the 
minimization of the quotient in equation is subject to only one orthogonality condition, 
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thus only one negative direction in the energy functional can be eliminated by satisfying 
this condition. 

The following simple strategy has been used to establish the stability. The eigenvalue 
problem was reformulated in the polar coordinates (p, 8) and the operators were expanded in 
Fourier series with respect to 6, which is done by the simple substitution: V 2 — > V 2 — n 2 /p 2 , 
where V 2 = d 2 + p~ l d p . Noticing that the orbital operators Ai >n = Ai(V 2 — > V 2 — n 2 / p 2 ) 
are ordered as follows Ai. n+ i > Ai j7l , we have checked for the negative eigenvalues of the 
first two orbital operators with n = 0, 1. If there are two or more negative eigenvalues (for 
instance, if A^o and Ai ; i both have one negative eigenvalue) then the solution is unstable; 
if there is only one such of A 10 then one can apply the VK criterion. From our numerical 
simulations it follows that Aj_ j. is always positive, while the operator A 10 has one negative 
eigenvalue or none (the latter case corresponds to the defocusing effective nonlinearity in 
the coupled-mode system with the transverse potential, see sections [IV Bl and W \ . 

When the operator Ai does not have negative eigenvalues at all the solution is uncon- 
ditionally stable. In this case, the formal threshold |p = of the stability criterion is the 
point where an additional zero mode of the operator Aj. ( q appears. Indeed, at the threshold 
point we have 

due to the identity 




which follows from differentiation of the stationary coupled-mode system, Aq(U, V) t = 0, 
with respect to p. Hence, in this case, on one side of the VK threshold the operator A 10 
is positive and the solution is unconditionally stable, while on the other side there is one 
negative eigenvalue and the VK criterion applies (this corresponds to a change of the effective 
interaction from the repulsive to attractive, see sections ll V Bl and IV jl . 

Note that the above approach allows us to decide on the stability of the stationary solu- 
tions to the coupled-mode system without numerical solution of the full eigenvalue problem 
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IV. SMALL-AMPLITUDE AND ASYMPTOTIC SOLUTIONS 



In this section we study two limiting cases of solutions of the coupled-mode system (|15|): 
the solution with vanishing amplitude (i.e., the bifurcation from zero) and the asymptotic 
solution with the chemical potential taking large negative values. There is an essential 
difference in the bifurcating small-amplitude solutions in the coupled-mode systems with 
and without the transverse potential, while the asymptotic solution, though also possessing 
some minor difference between the two cases, can be studied for both cases simultaneously. 
Accordingly, the section is divided into three subsections. First we study the soliton bifurca- 
tion from zero, i.e. the coupled-mode system without the transverse potential, in subsection 
II V Al In subsection II V Bl the same bifurcation is considered in the coupled-mode system 
with a general transverse potential. We derive the asymptotic solution in subsection IIV CI 
using, for simplicity, the parabolic transverse potential. 

A. The soliton bifurcation from zero 

Though the linear coupling of the repulsive and attractive NLS equations breaks the 
scale- invariance, the 2D-soliton solutions are always unstable. This conclusion follows from 
a comparison of the soliton bifurcation from zero, i.e. \x — > fiut (which is the subject of this 
section) and the asymptotic solution when /i — ► — oo (which is considered in section HV Cj) . 
Moreover it is supported by the direct numerical solution of section IV1 

The stationary coupled-mode system (when the transverse potential is flat) reads 



First of all, the most important information on the existence of solitons is provided by the 
dispersion law of the linearized system, which has two branches in the case of the coupled- 
mode system (J2EJ). Setting U = U e~ xp and V = V e~ xp , in the limit of vanishing U and V 
we obtain: 



fxU + V 2 U + U 3 + kV = 0, 
(p-e)V + V 2 V-clV 3 + kU = 0. 



(28b) 



(28a) 



A? 




(29) 



It is easy to see that the following inequalities hold: fi\ < < and /ii < e < /U 2 . 
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The same dispersion relations appear also in the ID well jl9|, their universality 

is due to the fact that the limiting values of the chemical potential are determined by the 
energy difference in the double-well trap and the trap asymmetry: 

= T~, o - o — I » ^2 = o + n • ( 30 ) 



l + x 2 2 J hw± ' \l + x 2 2J tko A 

In the ID case the positive solitons (UV > 0) bifurcate from zero at the lower energy level 
p = pi, while the non-positve ones (UV < 0) from the higher level p = p 2 - This property 
also holds in the 2D case. This can be shown as follows. First of all, only one branch of 
soliton solutions may correspond to each branch of the dispersion law. Second, there is 
a point po on the positive real line such that V 2 L r (p ) = 0. Indeed, setting £ = hip we 
have V 2 £/(p) = e~ 2 ^^^, whereas, for the solution with a finite Z 2 - n orm, the first derivative 
% = p^f- tends to zero as £ — >• ±oo (i.e., when p — >• or p — >• oo). Considering equation 
(|28aj) at p = po we obtain (pU + U 3 + kV)\ Po = 0, thus p < for UV > (recall that 
k > 0), i.e. the positive solitons belong to the Ai-branch of the dispersion law, while the 
non-positive ones to the A2-branch. 

Consider now the positive solitons with vanishing amplitude, i.e. in the limit p — * p^f = 
Pi. Set p = pi — e, with e — > and suppose that in this limit U = 0(e p ) and V = 0(e q ) 
with some p, q > 0. From system (|28|) we necessarily get p = q. Further steps are essentially 
the same as in the ID case First we expand V in the power series with respect to U 



and its derivatives, using for this goal equation (J2I 

V = k(3U - an 3 (3 4 U 3 + 3a 2 K 5 (3 7 U 5 + k(3 2 V 2 U + 0(e s ), (31) 

where (3 = (e — p) _1 and s = min{3p + l,p + 2, 7p}. Here /3 should also be expanded with 
respect to e. Second, the result is substituted into equation ()28a|) and we get 
K 2S \ ( - —,2jj , ( 1 clk 4 \ 3 3a 2 K 6 5 



- e \ 1 + T4) u+ \ 1 + A) vu + 1 1 " A ) u + ~W U = 0(e,) ' (32) 

From the coefficient at the cubic term and the definition of p 2 one immediately concludes 
that the soliton bifurcation from zero is possible if and only if e > e cr , where e cr is the same 
as in the ID case lfl| : 

£ci = K (a 1 / 4 - a" 1 / 4 ) . (33) 

For e > e CT one can drop the quintic term from equation (|32~j). thus the effective bifurcation 
equation is the canonical 2D NLS equation except for the coefficients. In this case the 
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condition that all terms have the same order requires that p = 1/2 and V 2 ~ e, i.e. there is 
a new length scale £ = \ft~p. The soliton solution reads 

U = v^4 {R^p) + tQ ^ tp) + ( e 2)} , Ao = (l + ^ ' (\ - ^ 7 . (34) 

Here is the Townes soliton. 

While the soliton amplitude approaches zero, the corresponding number of atoms has a 
finite, limit. Indeed, we have 

N u = j d 2 pU 2 = A 2 N th + 0(e), (35) 

with N t h being the threshold for collapse (|T9|) in a single 2D NLS equation. Thus, the 
bifurcation from zero is discontinuous in the 2D case. This is quite dissimilar to the ID 
case, where only on the boundary e = e cr the soliton bifurcation from zero is discontinuous 

Q. 

On the boundary e = e CT the effective equation is the quintic NLS equation and p — 1/4, 
exactly as in the ID case. However, in contrast to the ID case, the bifurcation from zero is 
singular. Now U = e l/i B (F (^/ep) + eF^y/ep) + C(e 2 )), with some B and F (£). Hence, 
the number of atoms N u oc e -1 / 2 as e — > 0. The bifurcating solitons are unstable in this case 
by the VK criterion. 

For e > e cr , to decide on the stability of the soliton solutions near the bifurcation point 
one has to find the derivative of the number of atoms with respect to e in the limit e — > 0. 
To this end the soliton solution up to the order e 3//2 is required. The necessary higher-order 
expansion for V in terms of U reads 

V = k(3U - aK 3 (3 4 U 3 + 3a 2 K 5 f3 7 U 5 + k(3 2 V 2 U + k(3 3 V 4 U - 2aK 3 (3 5 V 2 U 3 + 0{e 7/2 ). (36) 

Substituting this expression into equation ()28aj) one can get an equation for the correction 
G — G{£) to the soliton solution (J33j) . The result is 

&G = (1 - V 2 - - 3R 2 )Q = C 2 R, (37) 

with the operator £2 defined as follows 

H^r^-^+f^H^f^i < 38) 
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Now, we can find the total number of atoms N = N u + N v . We have 

N u = A 2 N th + 2eA 2 J d 2 £RG + C(e 2 ). (39) 
The v-component of the soliton up to the order e 3//2 follows from equations (|34jl and (|36j): 
V = V~eA - ! R + e (-(V}-l)R- A 2 ^R 3 +g]+ 0(e 2 ) ) . (40) 

V-1 I 4 ft / J 

Thus the number of atoms iV„ is given as 



2 2 

N v = A 2 ^N th + 2e^A 2 l + 0(e 2 ), (41) 

ft ft 



where we have denoted 



i s / d . £ - fl (l ( v|- - ^ + e) = / d'f _ 2 (i + rf 

Relation ()21|) for the Townes soliton has been used to simplify the above expression for X. 

The scalar product of R and Q, which enters the expression for the number of atoms, is, 
in fact, equal to zero. Indeed, using equation (pT7)) we get 

-i. 



d*£R8 = / d^R£i 1 £ 2 R 



but 



£r 1 « = ~(^ + fV r fl). (42) 

(Equation can be obtained by taking the derivative of the stationary 2D NLS equation with 
respect to chemical potential /i at // = — 1.) Therefore 

y d 2 f ir& = -i y d 2 f (RC 2 R + -iV^RCzR)^ = 

via integration by parts in the second term. Hence, the total number of atoms assumes the 
following form 



l2/\r A^A2 K f-> , a/?2 \1 \ v rM , .2, 



N — N u + N v — I 1 H — 2) ^o^th - 4eA 2 — 1 + — ^ iV th + 0(e 2 ). (43) 
V ft/ ft \ ft / 

Clearly, (for a > 0) in the vicinity of the bifurcation point /i — /ii we have |p = — ^ > 
which renders the two-dimensional bifurcating solitons unstable in contrast to the stability 
of the similar bifurcating solitons in one spatial dimension Q]. 
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B. Bifurcation from zero in the presence of transverse potential 



We have seen that the soliton bifurcation from zero is always discontinuous due to the 
fact that the solution has a new length scale £ = y/ep. If the transverse parabolic potential 
is taken into account (i.e., when it is not flat) the bifurcating solution has the length scale of 
order 1 (i.e., the order of the oscillator length). The stationary coupled- mode system with 
the parabolic transverse potential can be cast in a form analogous to that of system (|28j) : 

ujU + VU + U 3 + kV = 0, (44a) 
(uj-e)V + VV -aV 3 + kU = 0, (44b) 

where we have introduced the operator T> = V 2 + 2 — p 2 and a shifted chemical potential 
lj = p — 2. Note that T> < with T>e~ p2 l 2 = 0. System f)44j) in the limit of a vanishing 
amplitude solution, U = At~ p2 l 2 and V = Be~ p2 / 2 , with A, B — > 0, gives two boundary 
values of u which coincide with the limiting chemical potentials of section HV Al 2 = pi^- 
Moreover, the amplitudes are related as follows B X2 = — (^1,2/^)^1,2- Therefore, the positive 
solution AB > bifurcates at oj = Ui, exactly as in the soliton case. 

To study the bifurcation in detail let us set to = to\ — e with e — > (here e can be negative). 
Further steps in the derivation of the leading order equation for the bifurcating solution are 
formally the same as those in section HV Al one has only to substitute p —>■ oj and V 2 — > T>. 
For instance, we have U ~ V ~ |e| p , where p > 0. For V we obtain the expression formally 
equivalent to that of equation (J3T|) . 

V = Kj3U - aK 3 (3 4 U 3 + k(3 2 VU + 3aV/ff/ 5 + C(|e| s ), (45) 

with j3 = (e — oj)" 1 and s = min{3p + l,p + 2, 7p}, while U satisfies the equation 

1 + -2 V + f 1 + -2 ) VU + f 1 " V 3 + ^- U * = °^\ S )- (46) 

vi) \ viz \ ih ) i4 

In the derivation of equations (J43j) and ()46|) it is assumed that T> ~ e, an analog of V 2 ~ e 
of the previous section, though the reason is different: the operator T> is small because of its 
discrete spectrum and the fact that the solution bifurcates from the ground state (with zero 
eigenvalue). The critical zero-point energy difference e CT (133)) delineates the regions of the 
defocusing and focusing cases; though for e < e CT the effective equation (j4Tlj) is defocusing, 
it has a localized solution thanks to the external potential (in the operator T>) . 
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To make our approach general we will use only two properties of the operator D, namely 
that it is non-positive (with zero being an eigenvalue) and that it has a discrete spectrum. 
The case of e — e cr is a special case of the bifurcation from zero, exactly as for the soliton 
bifurcation of the previous section. Consider first e 7^ e CT , i.e. when the cubic term in 
equation (J46j) has a non-zero coefficient (positive or negative). Hence, p = 1/2 in this 
case. Defining a new dependent variable F by setting U = y/\e\\Ao\F(p), with Aq given by 
equation we obtain for F the equation 

-F + e^VF + aF 3 = 0(\e\ 3/2 ), a = sgn {e(e - £ cr )} , (47) 

which has a non-zero solution in the leading order e° only if a > 0, i.e. when the sign of e 
is equal to that of e — e CT . Equation (J47|) allows one to obtain two leading orders of F in 
the expansion with respect to e: F = Fq + eF% + 0(e 2 ). The simplest way to get them is to 
invert the operator 1 — e~ l T>: 



(l-e-^H^oi + e^^^ (48) 

71=1 

where we have used the eigenvalues, — A n (A n > for n > 1), and the eigenf unctions, \4> n ), 
of the operator T>: T>\cf) n ) = — \ n \4> n ). Now, multiplying equation (|4*7j) by the operator from 
equation (}4*%j) and collecting the successive orders of e we obtain 

F = (tir 112 U-Jh E **{P) + ^o)- V2 E { ^^Mp) + 0(e 2 ). (49) 

Here (...) denotes the integral J d 2 p(. . .). 

Let us find the number of atoms corresponding to the bifurcating solution. We need just 
the leading order, since, the bifurcation is continuous and the number of atoms tends to zero 
as uj — * uj\. We get 

N = J d 2 p(U 2 + V 2 ) = \e\(4)- 1 (l + £) Al + 0(e 2 ), (50) 

where we have used that in the leading order V = -j^U + 0(|e| 3 / 2 ). Now, for e > we have 
|p = — < 0, i.e., the solution is stable for e > e CT by the VK criterion (since there is 
only one negative eigenvalue of the operator Ai pijl ). In the case e < e CT the solution is 
unconditionally stable (the operator Ai is positive). Unconditional stability in the latter case 
can be explained by the fact that the effective equation pUj) is the defocusing NLS equation 
with the external potential, hence the ground state solution is unconditionally stable. 
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Consider now the special case of the bifurcation from zero, when e = e cv . The leading 
order nonlinearity is quintic in U and we get p = 1/4. Defining a new variable F by setting 
U = \e\ 1 ^ 4 B F(p), where Bq = ( 1 + ^ J 3 ^ 2 6 , we get the effective equation 



which allows one to obtain the leading order F = Fq + 0(e). It is clear that in this case 
e > 0. We get F Q = (^) _1/4 o (p)- Finally, using V = —U + C(e 3 / 4 ), we obtain 



Therefore in this case the solution is also stable by the VK criterion (there is only one 
negative eigenvalue of the operator Ai (J2U), similar as in the previous case). 

We see that the confining transverse potential allows for stable small-amplitude solutions. 
Moreover, for the zero-point energy difference below the critical value, e < e cr , the effective 
equation for the bifurcating solution is the defocusing NLS equation, whereas in the opposite 
case it is the focusing cubic or quintic NLS equation. 

C. Asymptotic solution of the coupled-mode system for fl — ► — OC 

In the previous two sections we have considered the bifurcation from zero. To complete the 
consideration, one has to study also the other limit of the chemical potential, i.e. /i — ► — oo. 
This can be done in a unified way for the coupled-mode system with or without the transverse 
potential. The reason is that in both cases there is the same new length scale and the solution 
is approximated by the Townes soliton. Therefore, we will use the stationary system (f4*4*j) 
to study the asymptotic solution. Thus we set u> = — e _1 with e — > 0, where uj = \i — 2 as 
in the previous section. Finally, we will restrict the consideration to the case a > 0, which 
is the most interesting one, and assume that the transverse potential is parabolic (which is 
not an essential requirement, but is convenient for calculations). 

Analysis of the system f)44j) reveals that in the limit u — > — oo the leading order of the 
solution is as follows: U = 0(e~ l l 2 ) and V = 0{e 1 ^ 2 ). Indeed, the leading order of V cannot 
be greater then that of U otherwise there is no localized solution (for a > 0) in the limit 
e — ► 0. Equation (|44ajl for U has a non-trivial solution only if T> — (9(e _1 ). The only way 
to satisfy the latter is to require that V 2 ~ e -1 , i.e. there is a new length scale £ = e _1 / 2 p 




F + e~ l VF + sgn(e)F 5 = 0(|e| 3/4 ) 



(51) 




(52) 
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(the effective length of the solution tends to zero as uo — > — oo and the external potential 
is negligible). We will need to compute U up to the order e 3//2 and V up to the order e 1 / 2 . 
Setting 

U = e- 1 /2( f/ (0)( e -l/2 p) + eE /(l)( e -V2 p ) + ,2^(2)^-1/2^ + 0(e 3 ))? 

V = e 1 /2 ( y(0) (e -i/2 p) + o(e))) 
and expanding system in the series with respect to e, we obtain U^ ' = -R(£) and 

AC/«(0 = 2it;(0, (53a) 

At/ (2) (0 = 2f/ (1) (0 - £ 2 J2(0 + 3i?(0t/ (1) (0 + (53b) 

(l-V£)V<°>(0 = «fl(0. (53c) 

where the operator £i is given in equation (|37jl . Since in the leading order we have the 
Townes soliton, e is positive, i.e. the solutions with u — > oo are impossible (though for 
e < e cr , according to the results of section IIVB1 the curve iV = N(u) corresponding the 
solution of the coupled-mode system with the transverse potential initially enters the right 
half-plane u > ui, it eventually turns left and approaches — oo, see figure El in the next 
section) . 

By setting k = in equation (|53bjl one must obtain the derivative |p corresponding to a 
single NLS equation (with or without the transverse potential) in the limit \i — > — oo. In the 
case without the transverse potential we know that this derivative is zero due to the critical 
scale invariance. In the case of a single 2D NLS equation with an external potential the 
derivative is negative and the number of atoms approaches the collapse threshold (|19|) 
from below. "Switching on" the quantum tunnelling changes this behavior: for k > K cr , with 
some n CI , the derivative assumes a positive value and the number of atoms approaches 
the collapse threshold iV t h from above. Let us find the critical tunnelling coefficient. We 
have: 

N u = J d 2 pU 2 = N th + e 2 J d 2 f {U^ 2 + 2RU^ + G(e 3 ) 
= N th -e 2 j d 2 i [i 2 R 2 ]+e 2 K 2 j d 2 f U^(l - V 2 -)" 1 ^ + 0(e 3 ), 
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N v = J d 2 pV 2 = eV J d 2 £ R(l - V|)~ 2 i? + 0(e 3 ). 

In the derivation of these formulae we have solved equation (|53bjl for and ()53cj) for 
by the inversion of the corresponding operators and integration by parts. We have also used 
that the scalar product J d 2 £ RU^ is equal to zero, which is an immediate consequence of 
the fact that equation (|53aj) for does not involve the external potential (equation (}4*2*j) 
also can be used to establish this fact directly). Using equation (|53a|) and integrating by 
parts to get rid of the term with in the last integral in the expression for N u , we get 
the following formula for the total number of atoms 

N = N u + N v = N th - t 2 J d 2 £t 2 R 2 + £ 2 k 2 J d 2 f R(l - Vfj-'R + (9(e 3 ). (54) 

For the coupled-mode system without the transverse potential, the first integral on the 
r.h.s. of equation (JH3|) is absent. The second integral is positive: 



h = J d 2 £R(l - Vf)^R w 7.41. (55) 

Hence, in this case, the number of atoms always approaches the threshold for collapse N t ^ 
of a single 2D NLS equation from above. 

Consider now the case of the parabolic transverse potential. In this case the second term 
on the r.h.s. of equation (|H%|) is negative and 



d z a z R A w 13.82. 

Therefore, the two terms of order e 2 on the r.h.s. of equation (|54j) compensate each other 
at some k = k cv and the derivative |p changes sign in the limit fi — > — oo. The derivative 
itself is, however, of the order — 1/u 3 = e 3 and the change of sign is not visible numerically. 

The important conclusion from the asymptotic solution of the coupled-mode system for 
H — ► — oo is that a new scale appears and the w-component of the solution is approximated 
by the Townes soliton, while the amplitude of the f-component tends to zero. Hence, there 
are always solutions which suffer from the collapse instability. They have the total number of 
atoms approaching A^ t h as \x — > — oo. The collapse instability also appears in this limit in the 
one-dimensional coupled-mode system pjj], where for a large number of atoms the solution 
with —>■ — oo is, in fact, the collapsing ground state. However, in the two-dimensional 
system with the transverse potential another stationary solution becomes the ground state 
(see the next section). 
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V. NUMERICAL SOLUTION OF THE COUPLED-MODE SYSTEM. 



From the analytical study of the bifurcation from zero and the asymptotic solution for 
jj, — > — oo we can conclude the following. First of all, in the case of the coupled-mode system 
without the transverse potential, the soliton solutions exist for p < /ibif with /i bif = pi 
from equation ()29|) . The expansions of the total number of atoms for the bifurcating and 
asymptotic solitons are as follows: 



2 \ 



N = 1 + — A 2 N th + 4(p-p hi{ )A 2 — 1 + —K N th + (D(n- Pmf, P -> /ibif, (56) 

V ft/ \ ft / 

iV = iV th + — A + O^" 3 ), /i^-oo, (57) 

where v4 is defined in equation (JMj) . -R(£) is the Townes soliton and I\ is given by formula 
(|55|). Hence, we have iV(/ibif) > A^h for all values of the system parameters (we consider 
a > 0). Moreover, in both limits the derivative |p is positive. This gives a clear indication 
that the total number of atoms is a monotonic increasing function of the chemical potential 
and the soliton solutions to the system of linearly coupled focusing and defocusing 2D NLS 
equations are unstable in the whole domain of their existence. This conclusion was verified 
numerically. In figure Q we show the number of atoms vs. the chemical potential for the 
2D-solitons (in all figures we use the "scaled number of atoms" N, related to the actual 
number of atoms by formula (jl7)l . and the dimensionless chemical potential which pertains 
to the coupled- mode system (JUJ)). 

For the coupled-mode system with the transverse potential the behavior of the number 
of atoms as a function of the chemical potential is much more interesting. In this case, 
the total number of atoms satisfies equations (J5Uj) . (J52)) and (J54|) . Generalizing these to an 
arbitrary confining potential, we have: 



N 



iM-MbifK^)- 1 (l + ^) l-^r +0(/x-/i bi f) 2 , e^Bcr, 

31 V2 (58) 
+ 0(n - ^bif), e = e crj 



0^-/^(06) -1/2 



^2 / 1 I K 



N = N th + ^ + i«t) + 0(fi- 3 ). (59) 

Here (■■■) = J d 2 p and (f>o(p) is the ground state wave- function of the linear operator — V 2 + 
Kxt(p)> where V ext (p) (an even function) is the confining potential; Ii is given by equation 
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(IHH1) . while J e xt is determined by the quadratic term in the Taylor expansion of the trap 
Kxt(p) about p = (for the parabolic trap V ex t = P 2 we have I ext ~ —13.82). 

As shown in the previous section, for \i — > — oo, the w-component of the solution tends 
to the Townes soliton, while the amplitude of the w-component tends to zero. Hence, the 
solution suffers from the collapse instability in this asymptotic limit. From equation (|59|) 
it is seen that the total number of atoms approaches N^. However, one cannot conclude 
that the collapsing solution is the (unstable) ground state neither that there are no stable 
solutions to the coupled-mode system, which have the total number of atoms greater than 
iV-th- In fact, such solutions do exist and give the ground state of the system, stable with 
respect to collapse. Figure 121 illustrates an example of the stable ground state solution with 
the total number of atoms larger than N t ^, there a = 0.005, n = 2, e = —3 and N = 21.5 
(it corresponds to a point with fi = —2.31 on the curve of figure 0] below). Indeed, equation 
(I58|) clearly states that the bifurcating solution is always stable (by the VK criterion for 
e > £ cr and unconditionally otherwise, according to the discussion in sections IlIII and TlV Bj) . 

To determine the ground state of the system we have calculated the energy of the solution 
numerically. In the coupled-mode approximation, the energy of the condensate in the double- 
well trap is given as 

E = fo^i. (Jd^l^l 4 )' 1 ^ 
167r|ai n ^| 

H = J dV||V M | 2 + |Vt;| 2 + p 2 (|M| 2 + |t;| 2 )+£|t;| 2 -/€(^ + ^*)-^ + a^|. (60) 

Here 7i is the Hamiltonian of the dimensionless coupled-mode system. 

The existence of stable stationary solutions with large total number of atoms depends on 
the system parameters, principally on the tunnelling coefficient k and the zero-point energy 
difference £. In the asymptotic limit of weakly coupled equations, i.e., when fi — > — oo, 
the ground state is unstable with respect to collapse, similar as in a single NLS equation 
with the external potential. For large values of n and large negative £ one can expect the 
appearance of a new ground state due to strong quantum tunnelling through the barrier and 
competition of the attraction in the w-condensate and the lower zero-point energy for the 
atoms in the v-condensate. The stable solutions with a large total number of atoms were 
indeed found, for instance, for a = 0, k = 10 and e = — 10; the corresponding dependence of 
the total number of atoms on the chemical potential is given in figure El We have checked 
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numerically (by a numerical analysis of the spectrum of Ai (|24jl from section that the 
VK criterion applies. Therefore, the stable solutions pertain to the region where ^ < 0, 
close to the point of the bifurcation from zero. 

Figure El corresponds to the special case of a = 0, i.e. when the f-condensate is non- 
interacting quantum gas. However, for a > 0, when the f-condensate is repulsive, the curve 
A" = N(fi) exhibits similar behavior, see figure 0] 

The part of the curve A" = N(fji) with ^ < 0, corresponding to the stable solutions with 
a large number of atoms, also minimizes the energy when there is a local maximum with 
the total number of atoms higher than A" t h- This is illustrated in figure El where we give 
the equation of state, 7i = 7i(N) (i.e., the energy vs. the total number of atoms in the 
dimensionless variables). The equation of state given in this figure is characteristic for the 
2D coupled-mode system for a large region of values of the system parameters, when there 
is a maximum of the total number of atoms higher than Ath- Otherwise, the energy takes 
positive values and the minimum (zero) corresponds to the collapsing ground state. 

Figures El and |U correspond to the case e > e cr . A new feature appears in the opposite 
case, i.e. when e < e cr . Indeed, in this case the effective equation (see equation (|47)j) from 
section lIVBj) is defocusing close to the bifurcation point. This fact changes drastically 
the dependence of the total number of atoms on chemical potential: the curve A" = N(fi) 
enters into the region of \i > /ibif, i.e. to the right of the bifurcation point, see figure El 
the inset. The protruding part of the curve is analogous to that in a single defocusing 2D 
NLS equation with an external potential. The operator Ai appearing in the linear stability 
analysis of section IIHI is positive definite there. Hence, the solutions corresponding to this 
part of the curve A" = N(fj) are unconditionally stable. 

Moreover, there is a turning point bifurcation, where ^ = oo (see the inset in figure El)- 
This bifurcation corresponds to the change of sign of the effective nonlinearity in the coupled- 
mode system from negative to positive as one moves upwards along the curve starting from 
the bifurcation point. Accordingly, the lowest positive eigenvalue of the operator Ai j0 passes 
through zero and becomes negative to the left of the turning point. As there are no other 
negative eigenvalues, the VK criterion applies to the left of the turning point bifurcation. 

The equation of state TL = TC(N) corresponding to figure El is similar to that illustrated 
in figure El Accordingly, the total number of atoms in the ground state exceeds A^h by an 
order of magnitude in this case. 
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Finally, figure [7| illustrates the fact that the maximum of the total number of atoms 
achievable by the stationary solution is determined by the tunnelling coefficient k and the 
zero-point energy difference e (compare to figure OJ). In this case e < e CT (the turning point 
bifurcation as well as the protruding part of the curve are also present , but not visible) . The 
stability is again determined by the VK criterion, except for the extremely narrow region 
before the turning point bifurcation (where the solution is unconditionally stable). In this 
case, the total number of atoms is always smaller than the threshold iV t h- 

As \i — > — oo the curves N = iV(/i) and N u>v = N u>v (n) in figures EE 13 and [7| are similar 
to that of the "solitonic" curve shown in figure H We have confirmed that the asymptotic 
solution for // — > — oo indeed approaches the Townes soliton in its -u-component, whereas 
the f-component tends to zero. Thus in the limit fi — > — oo the two condensates are weakly 
coupled, similar to the one-dimensional case 

The stable ground state in the 2D coupled-mode system, which corresponds to the part 
of the curve N = N(fi) immediately after the bifurcation point, appears due to breaking of 
the scale invariance by the transverse trap. This is reflected in the estimate of the solution 
width which is of the order of the oscillator length of the trap (see section ITVBj) . It is seen 
that the larger share of atoms is gathered in the repulsive f-condensate. Hence, such a 
state is similar to the unusual bright soliton of the one- dimensional coupled-mode system 
[l^ . However, there is an important difference between the one-dimensional and the two- 
dimensional cases: in one spatial dimension the ground state always corresponds to the 
weakly coupled condensates and suffers from the collapse instability for a large number of 
atoms. 



VI. CONCLUSION 



Nonlinearity of the Gross-Pitaevsky equation is due to the atomic interaction which is 
essential for understanding the properties of the condensate. Control over the nonlinearity in 
the Gross-Pitaevsky equation allows for coupling of two condensates with the nonlinearities 
of opposite signs (the scattering lengths). This can be realized using a double- well trap with 
far separated wells. In the one-dimensional case, there are stable bright solitons with almost 
all atoms gathered in the repulsive condensate In the two-dimensional case, on the 

other hand, the Townes-type solitons in the system are always unstable due to the fact that 
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mode system, 
3- 



the vanishing amplitude 2D-soliton solution has a finite ^-norm, i.e. the bifurcation from 
zero always corresponds to a discontinuity in the dependence of the number of atoms on the 
chemical potential. This is quite dissimilar to the one-dimensional couplec- 
where the bifurcation from zero is continuous except for the boundary case J 

With the use of a parabolic potential the spatial scale of the solution is fixed by the os- 
cillator length. This allows for the stable stationary solutions (though they are not solitons) 
with large ^-norms, which represent the ground state of the system. The ground state is 
secured from the collapse instability by an energy barrier. Interestingly, this ground state 
solution may have the /2-norm, i.e. the number of atoms of the condensate in a double-well 
trap, much higher than the collapse threshold in a single 2D NLS equation (in some cases, 
the total number of atoms exceeds the collapse threshold by an order of magnitude). This is 
a new phenomenon, which pertains only to the two-dimensional coupled-mode system, since 
in the one-dimensional case, for a large number of atoms, the ground state corresponds to 
weakly coupled condensates, has a large negative chemical potential and suffers from the 
collapse instability A more detailed study of the properties of the energy barrier for 
collapse is relegated to Ref. j^ . 
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FIG. 1: The total number of atoms (solid line) vs. the chemical potential corresponding to the 
2D-soliton solutions of the coupled-mode system. The dashed and dotted lines give the number 
of atoms in the u- and u-condensates, respectively. The number of atoms here is reduced by the 
factor A, as in equation (|17|). The chemical potential is given in the units of hui±/2. 
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FIG. 2: Ground state of the coupled-mode system (stable with respect to collapse). Here the 
parameters are a = 0.005, k = 2, e = —3, fi = —2.31 and N = 21.5. The order parameters U and 
V are given in the units of y/~K/£± and the radial length p in the units of £±. 
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FIG. 3: The total number of atoms vs. the chemical potential (the solid line) for a = 0. The 
dashed and dotted lines give the number of atoms in the u- and u-condensates, respectively. The 
axes units are as in figure ^ 
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FIG. 4: The total number of atoms vs. the chemical potential (the solid line) for a = 0.005. 
Here e > e cr . The dashed and dotted lines give the number of atoms in the u- and v-condensates, 
respectively. The axes units are as in figure ^ 
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FIG. 5: The energy of solutions vs. the total number of atoms, corresponding to figure |IJ The 
almost straight line contains the ground state of the system and corresponds to the part of the 
curve N = iV(/i) in figure QJ where ^ < 0. Here the number of atoms is reduced by the factor A 
from equation (|17|). whereas the energy is given in the units of huj±A/2. 
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FIG. 6: The total number of atoms vs. the chemical potential (the solid line). The dashed and 
dotted lines give the number of atoms in the u- and u-condensates, respectively. Here e < e cr . The 
inset shows a section of the figure about the bifurcation point. The axes units are as in figure ^ 
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FIG. 7: The total number of atoms vs. the chemical potential (the solid line) for small values of 
the tunnelling coefficient and energy difference. Here e < e cr . The dashed and dotted lines give 
the number of atoms in the u- and u-condensates, respectively. The axes units are as in figure Q 
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